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ABSTRACT 

Because massive, low-metallicity population III (PopIII) stars may produce very powerful 
long gamma-ray bursts (LGRBs), high-redshift GRB observations could probe the properties 
of the first stars. We analyze the correlation between early PopIII stars and LGRBs by using 
cosmological N-body/hydrodynamical simulations, which include detailed chemical evolu- 
tion, cooling, star formation, feedback effects and the transition between PopIII and more 
standard population I/II (PopII/I) stars. From the Swift observed rate of LGRBs, we esti- 
mate the fraction of black holes that will produce a GRB from PopII/I stars to be in the 
range 0.028 < Jgrb <0.140, depending on the assumed upper metallicity of the progeni- 
tor. Assuming that as of today no GRB event has been associated to a PopIII star, we esti- 
mate the upper limit for the fraction of LGRBs produced by PopIII stars to be in the range 
0.006< fcB.B <0.022. When we apply a detection threshold compatible with the BAT instru- 
ment, we find that the expected fraction of PopIII GRBs (GRB3) is ~ 10% of the full LGRB 
population at z >6, becoming as high has 40% at z >10. Finally, we study the properties of 
the galaxies hosting our sample of GRB3. We find that the average metallicity of the galaxies 
hosting a GRB3 is typically higher than the critical metallicity used to select the PopIII stars, 
due to the efficiency in polluting the gas above such low values. We also find that the highest 
probability of finding a GRB3 is within galaxies with a stellar mass < 10 7 M Q , independently 
from the redshift. 

Key words: gamma-rays: bursts - Population EH; method: numerical. 



1 INTRODUCTION 

Gamma-r ay bursts (GRBs) are the m ost energetic explosions in the 
Universe (Zhang & Meszaros 2004), offering exciting possibilities 
to study astrophysics in extreme conditions and up to very high 
redshi ft (z), as shown by the detection of GRB 090423 at z = 
8.2 ISalvaterra et alfeOOStlTanvir et al.ll2009h . According to recent 
studies and the so called "collapsar model", long GRBs (LGRBs), 
with durations longer than 2 sec, are linked with the final evolution- 
ary stages of massive stars. In particular, a Wolf-Rayet star can pro- 
duce a long GRB if its mass loss rate is "small", which is possible 
only when the stellar metallicity is very low. Indeed, when metal- 
licities are below ~ 0.1 — 0.3 Z@, the specific angular momentum 
of the progenitor allows the loss of the hydrogen envelope while 
prese rving the helium core dWooslev&Hegen l |2006t | Fryer et al.1 
1999). Observations of high- z LGRBs can provide unique informa- 
tion about the first stages of gala xy formation in the ea rly Universe 
and on the reionization process (M cOuin n et al. 2009). Mor eover, 
it has been recently propo sed dMeszaros & Reesll2010tlToma et al.l 
l201CUSuwa& Ioka 20LLb that also very massive, metal-free first 
stars, the so-called Population III (PopIII) stars, can produce a 
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LGRB. The rotation rate has as well important implications for the 
evolution of the star and the possible production of GRBs. Indeed 
the stars should retain sufficient spin to power the collapsar b urst 
teromm & Loeb!l2006l : IStacv et al.ll201 lklchiappini et al.)l201 ll) . 

Indeed, a recent study bv lSuwa&iokalJIoT J) shows that the 
jet can potentially break out the stellar surface even if the PopIII 
star has a massive hydrogen envelope, thanks to the final long- 
lasting accretion of the envelope itself onto the central black hole 
(BH). Since GRB luminosity and duration are found to be very sen- 
sitive to the core and envelope mass, they are probe of the fir st lumi- 
nous objects at the end of the dark ages. iToma et all d20 id) calcu- 
lated the afterglow spectra of such PopIII GRBs based on the stan- 
dard external shock model, showing that they could be detectable 
with the Swift BAT/XRT. Anyway, the rate of GRB from PopIII, 
even if measurable, is very small. Thus, observations of PopIII stars 
might be possible only with a large area and with the fields of view 
of future missions (as e.g. the Fermi LAT instruments). 

Thus, GRBs associated to PopII I stars should be observable 
out to redshift ~ 30 (e.g. Iciardi & Loebll2000l : IToma et al.ll201ol 
and references therein) and might repres ent the most promising 
way to directly detect the very first stars (Eamb & Reichartll2000l : 
iBromm & Loebll2002l : lsalvaterra et alfcoich . 
In this paper, we compare the redshift distribution of LGRBs from 
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PopII/I stars with the one expected by PopIII stars up to high red- 
shift , by using state-of-the-art N-body/hydrodynamical simulation s 
(see lMaioll2008l : iMaio et alj|201ol . 1201 ll ; iMaio & Iannuzz3l201lh . 
including detailed chemical evolution of the Universe, gas cool- 
ing, PopIII-to-PopII/I transition, and feedback effects. To select the 
candidates for LGRBs, we extract the information for the age and 
metallicity of newly formed stars, and adopt the collapsar model 
for both PopII/I and PopIII stars. We build three samples of possi- 
ble progenitors with different metallicity thresholds, and give esti- 
mates of the PopIII GRB rate measured by Swift. We additionally 
compute the contribution of the LGRBs from PopIII progenitors 
with respect to the total LGRBs rate. 

The paper is organized as follows: in Sec. [2] we present the N- 
body/hydrodynamical chemistry simulations used in this work; in 
Sec. [3] we describe the method for selecting different progenitors 
of LGRBs; we describe our results in Sec. [4] then, we describe the 
main properties of the galaxies hosting our LGRBs from PopIII in 
Sec|U finally, we discuss results and give our conclusions in Sec. 

m 



2 NUMERICAL SIMULATIONS 

We perform N-body/hydrodynamical simulations, including 
cosmological evolution of e~, H, H + , H~, He , He + , He ++ , H2 , 
H+ D, D+ , HD, HeH+ dYoshida et all 120031 ; IMaio et alj|2006l 
20071 120091) , PopIII and PopII/I star formatiori TTornatore et all 



Table 1. Simulation parameters. 



2007) according to corresponding initial mass fu nctions (IMF), gas 
cooling from resonant and fine- structure lin es jMaio et al.ll2007h 
and feedback effects dSpringel & Hernquisdl2003l) . The transition 
from the PopIII to the PopII/I regime is determined by the under- 
lying metallicity, Z, of the medium, compared to the critical value 
10~ 4 Zp) (see e.g. ISchneider et alj|2003l ; iBromm & Loebl 



l2003l:ISchneider et alfcOOd and references therein). If Z < Z cr it, 
then a Salpeter PopIII IMF is assumed with mass range between 
100 M Q and 500 M Q ; otherwise, a standard Salpeter IMF is 
adopted with mass range between 0.1 Mq and 100 Mq, and 
SNII range betwe en 8 M Q and 40 M (see IMaio et alj 1201 ll ; 
IBromm et al]|2009l) . 

The chemical model follows the detailed stellar evolution of 
each SPH particle. At every timestep, the abundances of the dif- 
ferent species are consistently derived, according to the lifetimes 
of the stars and the yields in the given mass rang^U- Metal mixing 
is mimicked by smoothing the metallicities over the SPH kernel. 
Pollution is driven by wind feedback, which ca uses metal spread - 
ing over scales of several kpc at each epoch (e.g. lMaio et al.l201lh . 
We warn the reader, though, that, due to poorly known quantities 
in stellar evolution modeling, different assumptions can determine 
variations in the PopIII-SFR esti mates of more than one order of 
magnitude, as already s hown in IMaio et al. I fcoid, I20T lj), and in 
M aio & Iannuzzil {201 ll) . More specifically, the largest effects are 
due to the PopIII IMF assumed in the simulations. Many studies 
suggest the PopIII IMF to be top-heavy, with stellar masses of the 



Lifetimes are from 



Maeder & Meynet 
ll995|) for SNII, 



Thielemann et al. 



Padovani & Matteucc I i ll 9931) (but see also, e.g. 



\\9Si ). The yields used are f r om [W oosley & Weaveil 
van den Hoek & G roenewegen 119971) for AGB stars, 



20031) for SNIa jHeger & Wooslevl fcpol for PISN. For 
a more extensive and detailed discussion on different IMFs and metal yi elds 
and their effect on the enrichm ent history, see also IMaio et alj feplOh (in 
particular Fig. 6 and Fig. 7), and Maio & Iannuzzi ( 201 1 )(e.g. Fig. 5). 



Box side 


Gas particle 


Dark matter particle 


Maximum physical 


[Mpc/A] 


mass [Mq/H] 


mass[Mo//i] 


softening [kpc//i] 


100 


3 X 10 s 


2 X 10 9 


7.50 


10 


3 x 10 5 


2 X 10 6 


0.50 


5 


4 x 10 4 


3 x 10 5 


0.25 



order of hundreds solar masses (e.g. lLarsonli998llAbel et al.l2002t 
lYoshidaetaIll2004l) . However, there is also evidence for a more 
standa rd, low-mass PopIII IMF, with masse s well below ~ 100 Mr 



(e.g. lYoshidi] 120061; lYoshida et alj 12007]; ICampbell & Lattanziol 
l2008t lSuda & Fuiimotdl2010l ; lGreif et aljhoilh . The main differ- 
ence between these two scenarios is the lifetime of the first SN, 
which is ~ 10 6 yr in the former case, and ~ 10 s yr in the latter, 
and this results in a transition from the PopIII to the Po pII/I regime 
which takes place earlier in t he top-heavy IMF scenario ( IMaio et al.1 
l20irj ; lMaio & Iarmuzzill2oTlh . On the other hand, uncertainties in 
the critical metallicity for the PopIII/PopII-I transition, Z cr u, are 
not expected to be crucial, since the rapid pollution process boosts 
metall icities up to ~ 1 0~ 3 ' Zp, in few 10 7 yr after the first explo- 
sions jMaio et al.ll20 id) . 

In the simulation used here, the cosmological field is sam- 
pled at redshift 2 = 100, with dark-matter and baryonic-matter 
species, according to the standard cosmological model, with total- 
matter density parameter at the present firj.m = 0.3, cosmological 
constant density parameter f2o,A = 0.7, baryonic-matter density 
parameter fio,b = 0.04, expansion rate in units of 100 km/s/Mpc 
h — 0.7, spectral normalization as = 0.9, and primordial spec- 
tral index n = 1. We consider three cubic volumes with comov- 
ing sides of 100 Mpc/ft, 10 Mpc//i, and 5 Mpc/ft, with correspond- 
ing gas mass resolution of ~ 3 x 10 s M Q /h, ~ 3 x 10 s M Q /h, 
~ 4 x 10 4 M //i, and dark-matter resolution of ~ 2 x 10 9 M /h, 
~ 2 x 10 6 M //i, ~ 3 x 10 5 M //i. 

The identification of the simulated objects (with their gaseous, dark 
and stellar components) was carried out by applying a friends-of- 
friends technique, with comoving linking length of 20% the mean 
inter-particle separation, and a minimum number of 32 particles. 
Substructures are identified by using a SubFind algorithm, which 
discriminates among bound and non-bound particles. The parame- 
ters used in the simulations are listed in Table [T] 

We show the analysis of the simulation outputs by reconstruct- 
ing the evolution of the star formation rates (SFRs) in our simula- 
tions and comparing with observed data in FigQ] We find that the 
total SFR of our larger simulation is in good agreement with the 
observational data, while the smaller boxes were not evolved be- 
low z ~ 5 because at these redshifts most of the power in structure 
formation is at larger scales. Nevertheless, we see that the total SFR 
converges at z ~ 5. In addition, the general trend is very similar, 
with the contribution from PopIII rapidly decreasing with redshift 
and reaching ~ 10 -4 the total SFR. Because of the higher res- 
olution obtained, differently from the 100 Mpc//i box, within the 
two smaller boxes, we marginally resolve the formation of miniha- 
los. Such primordial objects have typical masses below ~ 10 8 Mq 
and as small as ~ 10 s Mq and they witness the f orm ation of the 
first st ars. As detailed e.g. in lTornatore et al.ld2007l) and lMaio et al.1 
(2010), the first PopIII stars form in these small mass halos, of pris- 
tine composition, and their subsequent formation history depends 
on the details of metal enrichment and chemical feedback. For this 
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Figure 1. Star formation rate density, SFR [Mq yr _1 Mpc -3 ], as a func- 
tion of redshift for the simulation with comoving side box of 5 Mpc/fe 
(red-solid line), 10 Mpc//i (green-dashed line) and 100 Mpc/h (blue-dotted- 
dashes line) . For each of them, we show the total star formation rate density 
(upper lines) and the SFR of particles with Z < Zcrit (PopIII SFRs), i.e. 
that will produce a GRB3 (lower lines). Symbols with error bars are a com- 
pilation of observational data (Hopkins and Beacon 2006), the black-dotted 
line is the best fit for observational data (Li 2008b). 



reasons the smaller boxes are better suited to capture the star for- 
mation process at high redshift, when, according to the standard 
hierarchical scenario of structure formation, the galaxy population 
is dominated by objects with masses much smaller than the present 
ones. 



3 DIFFERENT PROGENITORS FOR GRBS 

We divide the possible progenitors in three subsamples: two asso- 
ciated with PopII stars (Z > Zcrit), with different upper cuts in 
metallicity, and one from PopIII stars (Z ^ Z cr i t ). As mentioned 
in Sec.[TJ recent studies on the final evolutionary stages of massive 
stars suggest that a Wolf-Rayet star can produce a LGRB if its mass 
loss rate is small. This is possible only if the metallicity of the star 
is low, and therefore the specific angular momentum of the progen- 
itor allows th e loss of the hydrogen envelope while pre serving the 
helium core (Woosley & Heger 2006; [Fryer et al. The loss 

of the envelope reduces the material that the jet needs to cross in 
order to escape, while the helium core should be massive enough 
to collapse and power a GRB. 

In order to select possible progenitors for LGRBs, we have ex- 
tracted from the simulations the information about the age and 
the metallicity of all the star particles. In particular, the metal- 
licity is computed as the ratio between the total mass of metals 
contained in star particles and the total stellar mass. According 
to the collapsar mode l, we consider young star particl es with age 
t < t c . = 5 x 10 7 yr dLapi et alj2008llNuza et al.l2007l) . We follow 
an approach similar to lCampisi et al.l d2009t) , and we will refer to 
our subsample in the following way: 

• GRB1, obtained by selecting star particles associated with 
PopII/I stars (Z > Z cr it), without upper metallicity cut; 



• GRB2, including star particles associated with PopII/I stars 
(Z > Z cr it), and with metallicity Z ^ 0.5Zq; 

• GRB3, defined by selecting star particles with metallicity Z ^ 

Zcrit- 

In particular, for the second subsample we adopt a metallic- 
ity cut of ~ 0.5 Zq because stronger metallicity cuts seem to be 
excluded by the a nalyses of the properties of observed GRB host 
galaxies at z < 1 (Mannucc i et alj2010l : ICampisi et alj|201ll) . 
In the following, we will briefly describe how we compute the rate 
of LGRBs (Sec.lOl. considering both PopII/I (Sec.[lIJ and PopIII 
(Sec. 13. 3b progenitors. 



3.1 LGRBs rate 

The formation rate of LGRBs in the simulations can be computed 
by: 



pGRB,i(z) = fcRB.i C,BH,i P*,i( z ) 



(1) 



where i indicates the 3 subsamples previously described, fGRB,i is 
the fraction of BHs that will produce a LGRB, £bh,% is the fraction 
of BHs formed per unit stellar mass, and p, t i(z) is the star forma- 
tion rate for the considered subsample, i. 

The observed rate (in units of [yr _1 sr ]) of LGRBs of the i-th 
subsample at redshift larger than z is then given by 



RGRB,i(> z) 



= lb 



dz 



, dV(z') 1 P _grbM) 
dz' 4tt 1 + z' 



dL'ipi(L') 



(2) 



where 75 is the beaming fraction, dV(z)/dz is the comoving vol- 
ume element, and tpi is the normalized LGRB luminosity function 
(LF). The factor (1 + z) _1 accounts for the cosmological time di- 
lation. The last integral gives the fraction of LGRBs with isotropic 
equivalent peak luminosities L > L t h, i.e. those LGRBs that can 
be actually detected by the satellite. The threshold luminosity Lth 
is obtained by imposing a 1-s peak photon flux limit of P — 0.4 ph 
s _1 cm -2 in the 15-150 keV band of Swift/BAT. 
There is a general consensus th at GRBs are jetted sources 
dWaxman etal.fl9 98; Rhoadl l 19971) implying fundamental correc- 
tions to their energy budget and GRB rates. We take into account 
this fact by considering that the observed LGRB rates should be 
corrected by the factor yi that is related to the jet opening angle 8 
by 7;, = (1 - cos 0) ~ 6 2 I2 dSari et al.lll998l) . Given the average 
value of 6 ~ 6° dGhirlanda et al.ll2007h , < -y b >~ 5.5 x 10" 3 , 
which is the value we adopt. 

Finally, in order to compute the observed rate we have to spec- 
ify the LF of LGRB, tpi(L), for the three LGRB subsamples con- 
sidered here. Although the number of LGRBs with good redshift 
determination has been largely increased by Swift, the sample is 
still too poor (and bias dominated) to allow a direct measure- 
ment of the LF. However, an estimate of the LGRB LF and of 
its evolution through cosmic times can be obtained by fitting the 
BATSE differential number count s and imposing the const r aint o n 
the bursts observed by Swift (e.g. Salyaterra & Chincarinil ( 2007); 
ISalvaterra et al.l d2009l):ICarnpisi et al.l d2010h . but see lButler et al.l 
(2010); IWanderman & Piranld2010l) for different approaches). 
In particular, for G RB 1-2 subsamples we adopt the LF obtained by 
ICampisi etaI1 d2010). who described the LF as 



t/,(L) oc ( — 



exp — — 



(3) 
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Table 2. Fraction of BHs that will produce GRBi, fGRBi ■ F° r * =1>2 two 
different values for the minimum stellar mass to die as a BH are considered. 
For i=3 two different values for the number of GRBs observed by Swift, 
NGRB Suli j t > are used to constrain Jgrbz ■ See text for details. 



fGRBi 


from BHs with 


from BHs with 


M > 20M Q 


M > 40M Q 


Igrbi 


~ 0.028 


~ 0.089 


Igrbi 


~ 0.044 


~ 0.140 




N GRB S 


wift 




500 (up) 


140 (up2) 


fGRB3 


~ 0.006 


~ 0.022 



where £ is the bright-end power index and L„ is a characteristic 
cut-off luminosity. In order to take into account possible evolution 
of the LGRB LF, the cut-off luminosity can be written as L* = 
Lo(l + z) a where Lq is the cut-off luminosity at z = 0. Since 
the SFR evolution of o ur GRB1 and GRB2 samples are similar 
to the Hostl sample in ICampisi et alj d2010h . we adopt here for 
both subsamples the following value for the LF free parameters: 
Lq = 0.3 x 10 50 erg s _1 , a = 2 and £ = -1.7, that provide a 
good description of the available data. For GRB3, we assume that 
the LGRB from PopIII stars are expected to be much brighter than 
PopII/ I LGRB. Since t ypical luminosities should exceed 10 53 6 erg 
s" 1 jToma et alj|20ich . we adopt L, = 10 54 erg s 1 constant in 
redshift and £ ~ —1.7. However, we consider a range of possible 
values with L, = 10 53 - 10 55 erg s" 1 and -1.5 < £ < -2.0. We 
will discuss the dependences of our results on these parameters in 
Sec.H 



3.2 LGRBs from PopII/I stars 

The LGRB formation rate (see definition in Eq. QJ can be ob- 
tained from the star formation rate in the i-th subsample once 
C,BH.i and faRB.i have been specified. The fraction of BHs 
formed per unit stellar mass can be computed by (bh,i = 
J^ 00 (^(m,)dro,/ Jqj m*(/>(m*)dm» where <j>(m*) is the 
Salpeter IMF and m m i n is the minimum stellar mass of stars dy- 
ing as BHs. Given the uncertainties on m m in we will consider here 
two cases with m m in = 20 Mq and m m i„ = 40 Mq, the latter 
consistent with the value assumed in the numerical simulations. 
Since not every BHs produces a LGRB, as described in eq.Q] we 
need to know the efficiency with which this happens, fGRB,i - This 
can be estimated by considering the observed rate of LGRBs de- 
tected by Swift. Indeed, as of today the Swift rate of LGRBs with 
redshift z > 1 is abouQ 57 yr _1 sr _1 . By imposing this value, 
we find f G RB,i =0.028 (0.089) and f G RB,i =0.044 (0.140), as- 
suming m m in = 20 Mq (m m i n — 40 Mq). These values corre- 
spond to having ~ 7 (~ 11) LGRBs every 1000 SNe for the GRB 1 
(GRB2) subsample. We point out here that to compute the /grb,i, 
we use the simulation with box sizes 100 Mpc/ft, since this is the 
one reaching lower redshifts. We then assume the same value of 
fGRB,i also for the other simulations. The adopted values are sum- 
marized in Table [2] 



2 Here, we assume that the fraction of LGRBs observ ed at redshift < 1 i s 
~ 20% of the whole sample of bursts, as estimated in Campi sTet aljfeoich . 



3.3 LGRBs from PopIII 

In a similar way, we have to specify the value of <^bh,3 and f G RB,3 
for GRB3. As already mentioned, the first stars have been assumed 
to be very massive, with masses > 100 Mq. Stars with masses 
between 14 — 260 Mq are expected to die as pair-instability 
supernovae (Heg er& Woosl ev 2002; Zeldovich & Nov ikov 1971, 
1999), leaving no compact remnant. 

Thus, only progenitors with masses in the range 100 — HOM© 
and 260 — 500M© will lead to the formation of a BH possibly 
triggering a LGRB event. For the assumed PopIII IMF we obtain 
Csh,3 -0.0032. 

The fraction of PopIII stars that will produce LGRBs is un- 
known. As of today, there is no strong evidence favoring the 
detection with Swift (or any other satellite) of a LGRB associated 
to a PopIII star. Indeed, even the most distant LGRB detected so 
far shares common prompt emission and afterglow properties with 
the lo w- and intermediate- z LGRB population dSalvaterra et al.l 
2009), suggesting that it likely originates from a normal PopII 
progenitor. Thus, we can set a firm upper limit on the fraction of 
PopIII stars triggering the LGRB event by imposing that no LGRB 
associated with a PopIII star has been detected by Swiyjfl that is: 

flGRB3(> 0) 1 (4) 

j Rgrb(>0) Nswift' 

where Nswift is the number of LGRBs observed by Swift. Thus 
an upper limit can be set by considering no PopIII LGRB in the 
whole sample of Swift bursts. At the time of writing this con- 
sists in 500, corresponding to /grb3 up = 0.006. A more gentle 
constraint can be obtain by considering that some PopIII LGRB 
can be hidden among those Swift bursts without redshift mea- 
surement. In this case, /grb3 up2 = 0.022 is obtained by im- 
posing no PopIII GRB present among the 140 LGRBs with mea- 
sured 2 (Table[2]>. These values correspond to have 1 GRB3 every 
< 2(7) x 10 5 M Q (Z < Z C rit) produced. We note that a part of 
the high-2 GRBs population might be missed by current follow up 
observations, as about 20% of all bursts appear as "dark", possibly 
hiding some high-z GRBs. However, there is increasing observa- 
tional evidence that most GRBs are "dark" because of dust extinc- 
tion dCenko et alj2009l ; lGreiner et al J201 ll) . 



4 RESULTS 

The redshift distributions of LGRBs of different subsamples ob- 
tained by integrating eq.|2]are shown in Fig. [2] Here we plot the 
evolution with redshift of observed cumulative rate for our subsam- 
ples of LGRBs in simulations with different box size. Dashed- red 
lines represent the rate of GRB1 (-Rgrbi), while dashed-dotted- 
blue lines the rate of GRB2 (_Rgrb2)- Note that the different choice 
of the minimum mass of the BHs does not alter the redshift evolu- 
tion of LGRBs. The solid-dark line is the upper limit for J?grb3 
case, obtained using /grb3 up2 , and dotted-line is the rate for the 
same sample computed using /grb3 up • We use the arrows to point 
out that the computed rate of GRB 3 are upper limits. In the bottom 
panel of each plot we report the ratio between the rate of GRB 3 
and the total rate, i.e. the sum of the rate of PopII/I GRBs and 
PopIII GRBs. The solid line refers to the case in which LGRB from 

3 Since the simulations stops at z = 1, when calculating i?GRB3(> 0) 
we assume that the SFR from PopIII stars is negligible at z < 1. 
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Figure 2. Top Panels: Observed rate for our samples of LGRBs in the three simulations with different side box. Dashed-red lines represent the rate of GRB 1 , 
dashed-dotted-blue lines the rate of GRB2, the solid-dark lines are the upper limits for the rate of GRB3 and dotted lines are the more restrict upper limit for 
the same sample (using /grb3 up antl /GRB3 up2 respectively). Bottom Panels: Evolution with redshift of the ratio between the expected rate of GRB3 and 
the total rate obtained summing the rate of the samples GRB3+GRB1 (upper solid-red line) and GRB3+GRB2 (lower dashed-dotted-blue line), for the three 
simulations with different size box. 



PopII/I stars comes from GRB 1 subsample whereas dot-dashed line 
to GRB2 subsample. Here we assume /grb3 up2 • 
If we consider the total rate of LGRB for the 100 Mpc//i simu- 
lation, i.e. including the contribution from both PopII/I and PopIII 
stars, the expected observed rate at redshift z >6 is about R to t (z > 



corresponding rate of GRB3 is below 0.017 (0 .06) sr 1 yr 1 



6) ~ 1 (fl/6 ) 2 /yr/sr. This is consistent with recent observational Trus is in agreement with a TOpll proge r 
constraints dPerlev et all 120091: iGreiner et all 1201 ll). as well as dSalvaterra et al.l2009l;rChandra et aTlhoiOl) 



for /g rb3 1it1 (/grb3 up2 ), and is consistent with ISalvaterra et all 
(2010) when the same value of /grb3 is adopted. At z > 6 
the expected fraction of PopIII LGRBs is ^ 10% and increases 
with redshift. However, at z ~ 8 it can not be larger than 20%. 
This is in agreement with a PopII proge nitor for GRB 090423 



with previous theoretical estimat es (Salvaterra & Chincarinill2007l ; 
ISalvaterraet"aTi r2009: Butl er et al .11201 Oh . 



At these redshifts, the 



To extend our results to very high z we turn now to our smaller box 
simulations. Indeed, as already discussed in Sec. 2, at these red- 
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shifts the cosmic SFR is dominated by low-mass objects that can 
be properly resolved only by our small-box simulations. We note 
that at z ~ 5 — 6 the results of the smaller simulations converge 
and are in good agreement with those obtained from the 100 Mpc/h 
side box run. As expected, the contribution of PopIII GRBs is more 
and more important as redshift increases, becoming dominant for 
z > 16. Thus, a LGRB detected at extremely high redshift would 
likely be originated from a PopIII star. 

These results are not very sensitive on the assumed LF for PopIII 
LGRB. In fact, we find that the rate of PopIII LGRBs at z > 6 
is within a factor of two if we vary the characteristic burst lumi- 
nosity in the range 10 53 — 10 55 erg s _1 and the slope in the range 
-1.5 < £ < -2.0. 



5 HOST GALAXIES OF GRB3 SAMPLE 

In Sec. [2] we describe how the host galaxies are obtained from the 
simulations. The number of expected LGRBs from PopIII in the k- 
th galaxy can be computed as N k , s = /grb,3 (bhs M,, 3 , fc (z), 
where Af»,3,k is the stellar mass in the fc-th galaxy that satisfies 
our selection criteria for the 3-rd subsample (see Sec. [3}. Here 
we are interested in characterizing the properties of PopIII GRB 
hosts at high z. To this end, we consider the results obtained in the 
10 Mpc/fa side box run (this choice is motivated by the fact that 
we need to resolve properly the PopIII hosts and to have sufficient 
statistics to compute average quantities). We then calculate the 
probability of having a GRB 3 in a host galaxy of a given stellar 
mass, M», and metallicity, Z, in different redshift bins. 
In Fig.[3]we show this probability, weighted for the total number of 
LGRBs, present in the full galaxy sample. Most of the simulated 
GRB3 hosts are expected to be found in a very well defined region 
of the M* — Z plane. These galaxies have always values of Z lower 
than the solar metallicity at all redshift, but typically higher than 
the critical value used to select PopIII stars (Zarit = 10 -4 Zq). 
This is because metal enrichment from massive primordial stars 
proceeds rapidly, and, as soon as they form within a galaxy, they 
quickly pollute the surrounding environment to values above 
Z crit (see e . g. IWise & Abell 120081; iGreif et alj|20ict IMaio et al.l 
l20ld 1201 ll; IMaio & Iannuzzil |201 lh . Despite this, pockets of 
metal free gas are fou nd within metal enriched galaxies (see e.g. 
iTornatore et al.l 120071) . where PopIII star formation can occurr. 
The details, though, depend on the underlying model of metal 
enrichment and on the level of turbulent mixing which cannot be 
properly resolved by cosmological simulations. 
GRB 3 also seem to reside typically in objects at the lower end of 
the stellar mass distribution, M» ~ 10 5 Mq, because these have a 
higher probability of hosting non polluted gas. On the other hand, 
they should rarely be found in objects with M* > 10 7 Mq. It 
should be noted that this result may depend on the resolution of 
the simulation. We expect in fact, the typical stellar mass of the 
host galaxy to decrease when minihalos are properly resolved. 



Given the expected relation between LGRB properties and the BH 
mass, PopIII LGRBs are predicted to be very brig ht and detectable 
with current and future instrument up to z > 20 1 Bromm & Loebl 



6 DISCUSSIONS AND CONCLUSIONS 

LGRBs are now detected up to extreme high redshift and are con- 
sidered promising tool to study the state of the Universe during and 
beyond the reionization epoch. Moreover, it has been recently ar- 
gued that even very ma ssive, metal-free first stars could in principl e 
trigger a LGRB event JWooslev & Hegedl2006t iFrver et al.ll 19991) ■ 



12004 iNaoz & Brombergl |2007t iToma et alj|2010h . If this is the 
case, the detection of PopIII LGRBs might represent the most 
promising way to dire ctly probe the nature of the very first stars 
isalvaterraetalfeOld) . 

In this work, we have used three N-body/hydrodynamical cosmo- 
logical simulations, with various box sizes (L=100 Mpc/h - 10 
Mpc/ft - 5Mpc//i), to constrain the number and redshift distribution 
of PopIII LGRBs and compare them to those of PopII/I LGRBs. 
To estimate the rate of LGRBs of various populations, we assume 
the collapsar model, and construct three samples of LGRBs: GRB 1 
are bursts produced by PopII/I stars, GRB2 produced by PopII/I 
stars with metallicity Z < 0.5Zq and GRB3 from PopIII stars 
(Z < 10~ 4 Zq). In order to compute the rate of PopII/I LGRBs 
at high z, we follow an observationally motivated approach, by 
adopting a LGRB luminosity function able to reproduce both the 
BATSE differential number c ounts and the redsh ift distribution 
of LGRBs detected by Swift dCampisi et al]l201(]h . We find that 
~ 1 yr _1 sr _1 PopII/I LGRBs are predicted to lie at z > 6, con- 
sistently with avail able observational constrains dPerlev et alj2009t 
Greineret al. 2011) as we l l as with most recent theoretical works 
(e.g.lSalvaterra et al.l2009tlBufler et al.l20ld:IWanderman & Piranl 
l2010l) . In order to constrain the contribution of PopIII LGRBs, we 
compute an upper limit on the fraction of PopIII stars that can 
power a LGRB, assuming that none of the bursts detected by Swift 
up-to-now comes from a PopIII progenitor. With this limit, we pre- 
dict that less than 0.06 sr" 1 yr _1 PopIII LGRBs should be detected 
at z > 6. PopIII LGRBs are found to be non-dominant at z < 10 
with respect to normal PopII/I LGRBs. However, as expected, their 
contribution increases with redshift and the probability to find a 
PopIII LGRBs at z > 8 is ~20%, reaching also ~ 40 - 50% at 
higher redshift, when the contribution of PopIII stars to the total 
star formation rate is maximum. This implies that if we observed a 
LGRB at redshift z > 10, it would have roughly the same proba- 
bility to be produced by a PopII/I or PopIII star. 
Taking advantage of the detailed description of the simulated galax- 
ies, we further investigate the properties of the possible hosts 
of PopIII-LGRBs at very-high redshift. We find that the average 
metallicity of the galaxies hosting a GRB3 is typically higher than 
the critical metallicity used to select the PopIII stars, due to the effi- 
ciency in polluting the gas above such low values. We also find that 
the highest probability of finding a GRB3 is within galaxies with a 
stellar mass < 10 Mq, independently from the redshift. 
They key assumption we have made throughout this paper is that 
the primordial PopIII IMF is top-heavy with mass range between 
IOOMq and 500Mq. However , there are many uncertainties about 
that and different studies (e.g . jYoshida 2006; Yoshid aeTaIll2007l: 
ICampbell & Lattanziol kood ; ISuda & Fuiimotol l2010h show that 
smaller-mass first stars could be born even at very early times. This 
would imply lon ger stellar lifetimes and two main pred ictable con- 
sequences (see lMaio et al.l2010l ; lMaio & Iannuzzil 20 llh : (i) a gen- 
eral delay in the cosmic enrichment history of roughly ~ 10 8 yr; 
(ii) a larger contribution to the PopIII star formation rate at high 
redshift. 

It is not the aim of the paper to investigate this alternative scenar- 
ios, but if this were true, as a consequence we would expect a cor- 
responding increment in the GRB3 rate at high z. 
Thanks to t he simultaneous observ ations from the Swift and Fermi 
satellites teissaldi. E. et al.l l201ll) . future detections of hieh-z 
GRBs might help discerning the contribution from PopIII stars, and 
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Figure 3. Probability of having a GRB3 in a host galaxy with given stellar mass and metallicity, in three different redshift bins. The contours refer to a 
probability of 100%, 75% and 25%, where the percentages refer to the contours from the outermost to the innermost, respectively. The horizontal line indicates 
the critical metallicity Z cr a . 



if at least 10 GRBs will be detected at z > 6 one should be a GRB3. 
At the present, though, the number of such high-z GRBs is too 
small and the available data not good enough to discriminate their 
origin. Nevertheless, the probability of detecting a high-z GRB3 is 
non-negligible, and it becomes comparable to the probability of de- 
tecting a GRB from PopII/I stars at z > 10. Thus, investigation of 
high-z GRBs with future missions (e.g. SWOM, EXIST, ORIGIN, 
Janus) could offer a realistic first opportunity to detect PopIII stars. 
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